I 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  PURI  All  01  STANDARDS  1RM  A 


title: 

author: 

date: 

sponsor: 

program  nos: 


NATURAL  FREQUENCIES  AND  MODE  SHAPES 
OF  CABLES  WITH  ATTACHED  MASSES 

S.  Sergev  and  W,  D.  Ivan* 

August  1980 

Naval  Facilities  Engineering  Command 

YF59.556.091t 01. 402 


s 


DTIC 

ELECT E 
°£C  1  ?  1980 


* 

California  Institute  of  Technology 

CIVIL  ENGINEERING  LABORATORY 

NAVAL  CONSTRUCTION  BATTALION  CENTER 
Port  Hueneme,  California  93043 


Approved  for  public  release;  distribution  unlimited. 


-  O 
JL  ^ 


Professor  Iwan  developed  the  mode  shape  solution  techniques 
described  here  under  contract  to  the  Civil  Engineering  Laboratory. 

The  authors  thank  California  Institute  of  Technology  graduate 
student  Shawn  Hall  for  his  efforts  in  preparing  the  apparatus  and 
conducting  the  experiments  to  determine  mode  shapes. 


Act.  J*>,  7t j 

;  hyCZJi-THl-ISfZ  i 

*■?..  .  j  „,,  I 

4— HW  *NHi'  •-  — ■  •  # 


Unclassified 


SECURITY  CLASSIFICATION  QT  THIS  FAOe  <Wh—\  On  Enfr»Q 


REPORT  DOCUMENTATION  PAGE 


MtoAf  nuOTI 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


TN-1583  i 


I*.  OOVT  ACCESSION  NO. 

DN787011 


%.  RECIPIENT'S  CATALOO  NIIMRER 

'-Aoi*  9(t  c 


4.  TITLE  (m*4  Sw4llit»J 


TURAL  JJIEQUENCIES  AND  MODE  SHAPES  | 
F  £ABLES  WITH  ATTACHED  JpASSES^  J 

1 


t.  Tm  OF  REPORT  A  PERIOD  COVERCO 

Not  final;  Jan  1978  —  Apr  1979 


t-  PERFORMING  O AO.  REPORT  HUNGER 


7.  AUTHOR 


••  CONTRACT  OR  GRANT  MONGERS  J 


io.  rnognam  element,  rroject,  tmk 

AHEA  0  WO  OK  UNIT  NUMBERS 


62759N; 


0.  RERTORMING  ORGANIZATION  NAME  AMO  ADDRESS 

CIVIL  ENGINEERING  LABORATORY 
Naval  Construction  Battalion  Center 
Port  Hueneme,  California  93043 


II.  CONTROLLING  OKFICE  NAME  AND  AOORESS 

Naval  Facilities  Engineering  Command 
Alexandria,  Virginia  22332 


IS.  HONITONINO  AOCNCV  NONE  0  AOORESVII  OlfMWiI  limm  CmnSIS  Olticm) 


01.402 


It.  SECURITY  CLOSE. 


Unclassified 

it*,  decloshkicotion/oownonooino 
SCHEDULE 


IS  OISTNIEUTION  STATEMENT  (ml  Rsrwmj 


Approved  for  public  release;  distribution  unlimited. 


IT.  OISTNIBUTION  STATEMENT  ( •>  «!•  aOalracI  AMwwW  In  Bl.c»  ».  II  MIMimiI  Aw*  RarwRJ 


19.  REV  VOROf  rCMftoiM  m  fulfil  i(R  If  Mrmarr  aid  lAAnrlfy  Ay  Mick  «•!•») 


Natural  frequencies,  mode  shapes,  strumming,  cables,  strumming  cables,  cable  dynamics, 
vibration. 


>0.  AMT R ACT  (Caatlmf  m  ravara*  ilR  It  nacaaiary  art*  IRndfr  Ay  Alack  waiAw) 

‘  An  algorithm  has  been  developed  to  calculate  mode  shapes  and  natural  frequencies  of 
taut  cables  with  attached  masses.  The  transcendental  equations  of  motion  are  solved  by  an 
iterative  technique  that  allows  accurate  calculation  of  extremely  high  mode  numbers.  The 
algorithm  has  been  implemented  as  a  FORTRAN  program  primarily  as  a  tool  in  determining 
drag  coefficients  of  submerged  strumming  cables;  however,  any  taut  cable  can  be  analyzed.  - 

(continued) 


oo,;: 


1473 


EDITION  Of  I  NOV  ts  IS  OBSOLETE 


Unclassified 

SECUBITT  CLASSIFICATION  OR  THIS  RAOE  <Bkww  Dm  CmmwS; 


Unclassified 

stcumry  CLAMiriCATtow  or  thu  c««  a«i«*t«j  _ 


20.  Continued 

>  To  assess  the  accuracy  of  the  program,  a  simple  experiment  was  conducted  to  determine 
the  natural  frequencies  and  mode  shapes  of  a  wire  with  attached  masses  driven  sinusoidally 
by  a  shaker.  The  algorithm  shows  close  agreement  with  the  experimental  dam. 


Library  Card 


Civil  Engineering  Laboratory 

NATURAL  FREQUENCIES  AND  MODE  SHAPES  OF  CABLES 
WITH  ATTACHED  MASSES,  by  S.  Seigev  and  W.  D.  lwan 
TN-1583  23  pp  ilhaa  August  19SO  Unclassified 


1.  Strumming  cabin 


2.  Drag  coefficients 


I.  YFS9. 3 36.091. 01. 402 


An  algorithm  has  been  developed  to  calculate  mode  shapes  and  natural  frequencies  of 
taut  cables  with  attached  masses.  The  transcendental  equations  of  motion  ate  solved  by  an 
iterative  technique  that  allows  accurate  calculation  of  extremely  high  mode  numbers.  The 
algorithm  has  been  implemented  as  a  FORTRAN  program  primarily  at  a  tool  in  determining 
drag  coefficients  of  submerged  strumming  cables;  however,  any  taut  cable  can  be  analyzed.  To 
aaaeaa  the  accuracy  of  the  program,  a  simple  experiment  was  conducted  to  determine  the 
natural  frequencies  and  mode  shapes  of  a  wire  with  attached  maaaes  driven  sinusoidally  by  a 
shaker.  The  algorithm  shows  cloae  agreement  with  the  experimental  data. 


Unclassified 


CONTENTS 


Page 


INTRODUCTION  .  I 

STRUMMING  OF  MULTISEGMENT  CABLE  SYSTEMS  .  2 

Cable  Dynamics  .  . 

Equations  of  Motion 
Boundary  Conditions 


Solution  Algorithm  for  Mode  Shapes  and  Frequencies 


Summary  of  Solution  Procedure 

EXPERIMENTAL  COMPARISON  .  8 

EXPERIMENTAL  VERSUS  COMPUTER  RESULTS  .  11 

CONCLUSION .  12 


looeaeionFor 


NTIS  GRMfcl 
UDC  TAB 

Unanr.ou.nced 

Justification 


By_ 

I  Avail  aud/ci- 
iDiat  t  special 


ft 


} 

i 


▼ 


9t«.UK) 


INTRODUCTION 


Many  ocean-founded  structures  rely  on  cables  for  support  or  anchor¬ 
ing.  Static  and  dynamic  modeling  of  the  interaction  of  a  structure  with 
the  sea  requires  calculation  of  the  loads  due  to  cable  drag.  Strumming 
(vortex-induced  vibration)  of  the  cables  increases  the  cable  drag  coef¬ 
ficient  and  causes  a  proportional  increase  in  the  total  drag  force. 
Several  methods  have  been  developed  to  predict  the  amplification  of  drag 
due  to  strumming.  These  methods  require  a  calculated  cable  mode  shape 
as  an  input  to  algorithms  that  determine  drag  amplification.  The  algor¬ 
ithms  have  been  developed  from  experimental  data  and  represent  the  state 
of  the  art  in  drag  amplification  prediction.  In  the  past,  matrix  methods 
have  been  used  to  determine  cable  mode  shape,  but  these  are  costly 
techniques,  subject  to  inaccuracy,  particularly  for  higher  modes. 

The  Civil  Engineering  Laboratory  (CEL),  under  sponsorship  of  the 
Naval  Facilities  Engineering  Command,  has  been  developing  computer 
models  for  analysing  the  response  of  cable  structures.  As  part  of  this 
project,  CEL  has  developed  an  iterative  solution  of  the  transcendental 
equations  describing  cable  mode  shape.  The  iterative  approach  is  fast, 
accurate,  and  can  easily  accommodate  a  variety  of  system  configurations, 
including  bodies  attached  to  the  cable. 

One  of  the  initial  cases  analyzed  by  the  new  iterative  technique 
waa  a  long  cable  with  380  attached  bodies.  This  cable  is  representative 
of  those  encountered  in  ocean  engineering  applications.  Such  cables  are 
often  excited  by  ocean  currents  to  vibrate  at  very  high  mode  nua&ers. 

The  calculated  mode  shape  for  mode  162  is  shown  in  Figure  1.  It  is 
apparent  that  the  mode  shape  is  complex  and  much  beyond  the  range  of 
intuition. 

Since  tne  mode  shape  calculation  has  a  significant  influence  on  the 
drag  calculation,  it  is  considered  essential  that  the  determination 
scheme  for  the  mode  shape  be  verified  by  comparison  to  experimental 
data.  Therefore,  an  experiment  was  carried  out  to  measure  natural 
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frequencies  and  mode  shapes.  The  description  of  the  iterative  algorithm, 
the  Measured  results  of  the  experiment ,  and  a  comparison  between  the 
experimental  data  and  the  computed  results  are  presented. 
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Figure  1.  High  mode  number  for  a  real  cable  system. 


STRUMMING  OF  MULTISEGMENT  CABLE  SYSTEMS 

The  amplitude  of  oscillation  and  the  effective  drag  coefficient  of 
a  strumming  cable  system  are  determined  in  the  following  analysis.  The 
natural  frequencies  and  mode  shapes  of  cable  oscillation  are  obtained  by 
an  iterative  substitution  algorithm  that  finds  the  solution  satisfying 
the  imposed  boundary  conditions  of  the  problem. 


Cable  Dynamics 

The  cable  system  considered  is  shown  in  Figure  2.  It  consists  of 
n  cable  segments  attached  to  n-1  masses.  The  cable  segments  have  an 
effective  mass  per  unit  length,  p^,  and  a  tension,  T^,  which  is  assumed 
to  be  constant  over  the  length,  JL,  of  the  segment.  The  effective 
attached  mass  (including  added  mass)  is  M^. 
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Y^x)  =  sin  a.  w  x  +  B.  cos  flr  ui  x  0  £  x  i  JL  (3) 

x  “  1 1 • •  •  ,n 

where  a±  =  yj  pi/Ti  (4) 

Boundary  Conditions 

In  addition  to  satisfying  the  equation  of  sot ion,  the  deflection  of 
each  cable  segment  aust  satisfy  certain  boundary  conditions.  These 
conditions  result  froa  the  geoaetric  conditions  iaposed  on  the  ends  of 
the  cable  assembly,  the  continuity  of  displaceaent  at  each  attached  aass 
and  the  balance  of  forces  at  each  attached  aass. 

Continuity  of  Displaceaent.  Displaceaent  aust  be  continuous  at 
each  attached  aass.  Therefore, 

Yi+1(°)  =  Y.(£.)  (5) 


Substituting  froa  Equation  3  yields 

B. =  A.  sin  a.  ui  £.  +  B.  cos  cr.u>  £. 
l+l  i  ill  ii 


(6) 


If  Ai  and  B^  are  known,  Equation  6  can  be  used  to  find  Bi+J. 

Force  Balance.  Figure  4  shows  the  forces  acting  at  each  attached 
aass.  The  balance  of  forces  at  each  of  these  points  gives 


T 

l+l 


!Zi+i 
9x 


x=0 


T. 

i 


dx 


=  M. 


x=£. 

i 


8tJ 


x=0 


(7) 
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Figure  4.  Force  balance. 


Substituting  from  Equation  3  yields 


T. a. A. -  (a.  A.  cos  a.  u>  i. 
l+l  i+I  i+I  vii  i  i 


-  a.  B.  sin  a.  u>  £.)  T.  =  -  M.  u>  B... 
li  i  i'  i  i  i+I 


i  ~  1 j  *  • « ^n 


Using  Equation  6  and  solving  Equation  8  for  Ai+1  gives 


A. 


1+1  «t+i  Ti+i 


(cl  cos  on  u>  -  IT  ID  sin  cl  ui  JL) 

(cl  sin  cl  w  *»  cos  *#  ij)  Bi|  (9) 


If  A^  and  are  known,  Equation  9  can  be  used  to  find  A^+j . 

Geometric  Boundary  Conditions.  At  the  left-hand  end  of  the  cable 
assembly,  the  displacement  is  assumed  to  be  zero.  Thus, 


Yj(0)  =  0 

Substituting  from  Equation  3  thus  implies  that 


(10) 


(ID 
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Since  the  scale  of  the  deflected  shape  of  the  cable  (node  shape)  is,  at 
this  point,  arbitrary,  let 


A1  =  1 


(12) 


With  Conditions  11  and  12  on  Aj  and  and  Equations  6  and  9  for  A^+j 
and  ,  all  subsequent  A's  and  B's  can  be  determined  provided  u>  (the 
natural  frequency)  is  known. 

Hie  system  must  satisfy  one  additional  boundary  condition  at  the 
right-hand  end  of  the  assembly,  where  the  displacement  is  again  assumed 
to  be  zero.  This  gives 


W  *  °  * 


w« 


(13) 


Equation  13  in  turn  implies  that 


Bn+1  *  0 


(14) 


The  values  of  ui  that  give  solutions  satisfying  Condition  14  are  the 
natural  frequencies  of  the  problem. 

Solution  Algorithm  for  Mode  Shapes  and  Frequencies 

If  tu  is  varied  from  zero  to  some  large  value  and  the  corresponding 

values  of  BQ+1  are  calculated,  the  result  will  be  as  shown  in  Figure  5. 

Each  point  for  which  Bq+j  =  0  represents  a  valid  solution  of  the  free 

oscillation  problem.  The  ui^  so  obtained  are  the  natural  frequencies  of 

the  system.  There  are  an  infinite  number  of  such  frequencies. 

The  mode  shape  associated  with  each  natural  frequency  w.  will  be 
Tkl  * 

denoted  by  Y^  (x).  Then, 


Y<k)(x)  = 


i(k) 


sin  a 


i“k 


x  +  B 


00 


COS  OL  X 


(15) 


6 


Figure  5.  System  natural  frequencies. 


Let  u»k  be  the  natural  frequency.  Then  the  deflected  shape  of 
the  cable  system  will  be  such  that  the  number  of  internal  zero  crossings 
(nodes)  is  equal  to  k  -  1.  The  node  number  of  a  particular  mode  shape 
can  therefore  be  determined  by  counting  the  number  of  internal  zeros 
associated  with  the  function  Y^j  i  =  l,...,n. 

Summary  of  Solution  Procedure 

The  solution  process  for  the  mode  shapes  and  frequencies  is  sum¬ 
marized  as  follows: 

1.  Assume  a  value  for 

2.  Let  Bj  =  0,  Aj  =  1. 

3.  Solve  for  B2,A2;  B^A^;...;  Bn,An;  Bq+1  from  Equations  6  and  9. 

4.  Check  for  Bq+1  =  0.  If  Bq+j  t  0,  compare  with  previous 
value  and  estimate  a  new  trial  value  for  u). 

5.  Go  to  step  2  and  repeat  until  Bq+j  is  less  than  some  prescribed 
value  or  the  change  in  ut^  is  less  than  some  prescribed  limit. 

6.  Determine  the  mode  number  by  calculating  the  number  of  internal 
zeros  of  the  mode  shape  Y^(x);  i  =  l,...,n. 
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EXPERIMENTAL  COMPARISON 


la  order  to  verily  the  Accuracy  of  the  algorithm,  a  simple,  easily 
conducted  experiment  was  designed.  A  taut  piano  vire  was  fixed  at  one 
end  and  excited  sinusoidally  at  the  other  end  by  a  shaker;  natural 
frequencies  were  determined  by  varying  the  excitation  frequency.  The 
experimental  configuration  is  shown  in  Figure  6;  wire  properties  are 
given  in  the  figure.  Two  cases  of  attached  bodies  were  considered: 

Case  1:  Four  equally  spaced  bodies,  each  with  a  mass  of 
18.3  x  10“5  slug 

Case  2:  Six  unequally  spaced  bodies,  each  with  a  mass  of 
12.9  x  10-5  slug 

The  objective  of  the  experiment  was  to  find  the  first  8  to  10 
natural  frequencies  and  mode  shapes  for  each  case  for  comparison  with 
calculations. 

Tension  was  measured  by  shaking  the  unloaded  wire  (no  attached 
masses)  in  the  5^  mode  (n  =  5).  Period  (sec)  was  measured  by  a 
digital  counter  accurate  to  0.01  ms.  Tension  T  (lbf)  was  then  calcu¬ 
lated  by  the  natural  frequency  formula  for  an  ideal  cable: 


Split-shot  lead  weights  (the  type  used  in  fishing)  were  attached  at 
appropriate  locations  with  a  small  amount  of  contact  cement.*  When  ten 
samples  of  shot  were  weighed,  the  mass  variation  (standard  deviation) 
for  a  given  size  of  shot  was  found  to  be  about  ±1%  of  the  mean  value 
cited.  The  size  of  the  masses  was  approximately  0.32  inch  in  the  x- 
direction. 

*The  cement  was  necessary  to  prevent  the  masses  from  shaking  loose 
during  the  tests;  its  weight  is  assumed  negligible. 
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Locations  x  measured  from  the  right. 

Ao  iB  Peak-to-peak  shaker  amplitude;  excitation  is  perpendicular  to  wire 


Except  for  some  very  high  frequency  noise  caused  by  air  rushing 
through  the  shaker's  air  bearing,  excitation  was  harmonic,  as  measured 
by  an  accelerometer  mounted  directly  on  the  shaker.  For  each  experimental 
run,  the  shaker's  displacement  amplitude,  Aq,  was  reckoned  from 

Displacement  Amplitude  =  - - — r  [Acceleration  Amplitude] 

(2  n  f)2 

The  antinode  amplitudes  recorded  in  the  data  have  not  been  corrected  for 
the  shaker  amplitude  (the  correction  is  usually  negligible).  However, 
if  the  correction  is  desired,  the  following  formula  should  be  used: 

A  ~A  a  (  X*nt  inode 

corrected  "measured  o  l  l 

The  shaker  frequency,  f,  was  monitored  by  a  digital  counter  capable 
of  measuring  the  period  to  within  ±0.01  ms.  Experimentally,  the  shaker 
frequency,  f,  was  decreased  from  some  arbitrary  frequency,  f#,  lying 
well  above  fQ,  to  f^,  at  which  point  a  sudden,  highly  visible  increase 
in  amplitude  occurred.  The  frequency,  f^,  was  then  recorded  as  the 
"natural  frequency." 

The  location  of  nodes  was  measured  by  eye  or,  in  the  rare  case  of 
very  small  amplitudes,  with  the  help  of  a  small  strip  of  paper  against 
the  vibrating  wire.  Obviously,  the  uncertainty  in  measurement  is  least 
where  the  nodes  are  well  defined;  i.e.,  when  the  surrounding  antinodes 
have  large  amplitude  and/or  are  close  together.  Thus,  depending  on  the 
mode  shape,  accuracy  varied  from  about  ±0.1  inch  to  about  ±0.75  inch. 

In  general,  of  all  the  quantities  measured,  node  location  was  the  most 
accurate . 

The  location  of  the  antinode  was  also  measured  by  eye.  Judgment  of 
antinode  location  by  this  method  is  very  difficult,  inasmuch  as  the 
antinodes  are  not  sharply  defined.  Every  effort  was  made  to  make  the 
measurement  as  accurate  as  possible  without  resorting  to  sophisticated 
equipment.  A  rule  was  placed  a  short  distance  below  the  wire  and  observed 
from  a  fixed  distance  above  the  wire.  A  strobe  light  running  slightly 


slower  or  faster  than  the  vibration  frequency  was  used  to  illuainate  the 
vibrating  wire;  this  made  the  wire  appear  to  oscillate  slowly  in  its 
mode  shape.  In  this  way,  the  amplitude,  A-e#gure(j,  which  is  reported  in 
the  data  tables,  was  read  to  an  estimated  accuracy  of  0.01  inch. 

Three  runs  of  mode  5,  Case  1  were  made  in  succession.  Two  runs 
were  as  identical  as  possible,  with  four  masses  equally  spaced.  For  the 
third,  the  mass  at  x  =  138.0  inches  was  moved  to  x  =  138.1  inches.  On 
the  basis  of  these  runs,  it  was  concluded  that  the  repeatability  of  the 
antinode  amplitude  measurements  was  of  the  order  of  5%,  and  could  be  as 
large  as  10%  under  circumstances  less  controlled  than  those  in  the  tests 
described.  It  was  also  observed  that  the  amplitude  ratios  were  fairly 
critically  dependent  upon  the  exact  location  of  the  masses. 

To  judge  the  bias  in  the  measurements  caused  by  the  asymmetrical 
excitation  (shaker  at  one  end),  several  cases  were  rerun  with  the  con¬ 
figuration  of  the  masses  mirrored  about  x  =  i/2.  The  predominant  effect 
of  asymmetrical  excitation  is  a  preference  for  certain  modes  over  others. 
For  example ,  in  one  run,  mode  10  was  easily  excitable,  but  not  excitable 
at  all  in  the  mirrored  run.  The  reverse  was  true  for  mode  8. 


EXPERIMENTAL  VERSUS  COMPUTER  RESULTS 

From  Tables  1  and  2  it  is  observed  that  the  agreement  between  the 
calculated  and  experimentally  determined  natural  frequencies  is  well 
within  the  range  of  the  expected  experimental  error.  The  node  and 
antinode  locations  also  show  good  agreement,  especially  for  the  lower 
mode  numbers.  Modes  of  order  eight  and  higher  were  rather  difficult  to 
observe  accurately  with  the  experimental  setup.  However,  the  results 
for  the  higher  modes  still  appear  to  verify  the  computer  model. 

Response  amplitude  was  compared  on  a  normalized  basis.  To  compare 
on  an  absolute  basis  would  have  required  that  both  internal  and  external 
damping  be  known  and  modeled  and  that  the  precise  nature  of  the  excitation 
be  known.  Since  this  infonation  was  not  readily  available,  the  nor¬ 
malized  response  was  used  for  comparison.  As  shows  in  the  tables,  the 
normalized  amplitudes  of  response  compare  favorably. 
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The  nature  of  the  node  shapes  of  the  first  15  nodes  for  the  six 
■ass  cases  is  indicated  in  Figure  7.  Modes  1  through  10  are  those 
presented  in  Table  2.  The  additional  nodes  shown  illustrate  that  the 
non-unifomly  distributed  nesses  cause  the  node  shapes  to  becone 
increasingly  conplex.  It  is  interesting  to  note  that  the  sinusoidal 
portions  of  the  node  shapes,  shown  clearly  in  node  15,  but  also  evident 
in  other  nodes,  represent  the  free  response  of  the  wire  unaffected  by 
the  attached  nasses.  This  figure  was  generated  by  the  conputer  algorithm. 


CONCLUSION 

A  new  conputer  algorithm  for  computing  the  natural  frequencies  and 
■ode  shapes  of  cable  systens  has  been  presented.  This  algorithm  is 
computationally  efficient  and  shows  excellent  convergence,  even  for  very 
high  node  numbers.  The  algorithm  is  also  quite  flexible  since  it  has 
the  capacity  for  treating  a  wide  range  of  different  systen  configurations. 

A  comparison  of  the  results  of  the  conputer  nodel  with  the  results 
of  an  experimental  study  of  a  vibrating  cable  systen  in  air  has  shown 
that  the  nodel  accurately  predicts  the  natural  frequencies,  node  and 
antinode  locations,  and  relative  response  amplitudes  for  all  nodes 
obtainable  experimentally.  The  extrenely  close  agreement  obtained  for 
■odes  up  to  nunber  10  does  not  insure  that  higher  node  shapes  are  correctly 
calculated,  but  it  does  give  increased  confidence  that  the  solution 
technique  can  be  validly  extrapolated  as  required  by  the  users  needs. 


Figure  7.  Continued 
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Computed  values. 

'Could  not  be  determined  accurately. 
Hot  excitable  by  experimental  setup. 
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SOUTHWEST  RSCH  INST  R.  DeHart,  San  Antonio  TX 
STANFORD  UNIVERSITY  Engr  Ub,  Stanford  CA 
STATE  UN1V.  OF  NEW  YORK  Buffalo,  NY 

UNIVERSITY  OF  HAWAII  HONOLULU,  HI  (SCIENCE  AND  TECH.  DIV.);  Ocean  Engmg  Dept 
UNIVERSITY  OF  ILLINOIS  URBANA,  IL  (LIBRARY) 

UNIVERSITY  OF  MASSACHUSETTS  (Heronemus),  Amherst  MA  CE  Dept 
UNIVERSITY  OF  NEW  HAMPSHIRE  (Corell)  Durham,  NH 
UNIVERSITY  OF  NOTRE  DAME  Katona,  Notre  Dame,  IN 
UNIVERSITY  OF  RHODE  ISLAND  Narragansett  RI  (Pell  Marine  Sci.  Ub.) 

UNIVERSITY  OF  SO.  CALIFORNIA  Univ  So.  Calif 
UNIVERSITY  OF  TEXAS  Inst.  Marine  Sci  (Library),  Port  Arkansas  TX 
AMETEK  Offshore  Res.  &  Engr  Div 

COLUMBIA  GULF  TRANSMISSION  CO.  HOUSTON,  TX  (ENG.  LIB.) 

EG&G  WASH  ANALYTICAL  SERV  CTR,  INC  Rockville,  MD 
EXXON  PRODUCTION  RESEARCH  CO  Houston  TX  (A.  Butler  Jr) 

GOULD  INC.  Shady  Side  MD  (Ches.  Inst.  Div.,  W.  Paul) 

GRUMMAN  AEROSPACE  CORP.  Bethpage  NY  (Tech.  Info.  Ctr) 

HONEYWELL,  INC.  Minneapolis  MN  (Residential  Engr  Lib.) 

MOBIL  PIPE  LINE  CO.  DALLAS,  TX  MGR  OF  ENGR  (NOACK) 

OCEAN  ELECTRONICS  APPLICATION,  INC.  (Softley)  Key  Biscayne.  FL 
SANDIA  LABORATORIES  Library  Div.,  Livermore  CA 
TRW  SYSTEMS  CLEVELAND,  OH  (ENG.  LIB.) 

WESTINGHOUSE  ELECTRIC  CORP.  Annapolis  MD  (Oceanic  Div  Lib,  Bryan) 

WESTINTRUCORP  Egerton.  Oxnard.  CA 

WM  CLAPP  LABS  -  BATTELLE  DUXBURY,  MA  (LIBRARY) 
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